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Abstract 

We study the interplay between topological observables and chiral and Higgs tran- 
sitions in lattice scalar QED with quenched fermions. Emphasis is put on the chiral 
transition line and magnetic monopole percolation at strong gauge coupling. We 
confirm that at infinite gauge coupling the chiral transition is described by mean 
field exponents. We find a rich and complicated behaviour at the endpoint of the 
Higgs transition line which hampers a satisfactory analysis of the chiral transition. 
We study in detail an intermediate coupling, where the data are consistent both with 
a trivial chiral transition clearly separated from monopole percolation and with a 
chiral transition coincident with monopole percolation, and characterized by the 
same critical exponent v ~ 0.65. We discuss the relevance (or lack thereof) of these 
quenched results to our understanding of the xU(f>4 model. We comment on the 
interplay of magnetic monopoles and fermion dynamics in more general contexts. 
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1 Introduction 



Lattice gauge models with charged scalars and fermions have been studied 
both analytically and numerically for some time [1,2]. Their phase diagrams 
display a rich structure with chiral and Higgs transitions, as well as monopole 
condensation and percolation, and other topological excitations. 

In this paper we study the interplay of chiral symmetry breaking and monopole 
percolation in scalar QED with quenched fermions. Abelian models can be 
rewritten by a duality transformation as theories of topological excitations, 
which help in determining their phase structure and critical behaviour [3]. A 
classic example is the confinement transition of the pure compact U{1) lattice 
model, where monopole condensation was identified long ago as the underlying 
physical mechanism [4]. It is then interesting to investigate whether magnetic 
monopoles affect fermionic theories, and, in particular, the physics of chiral 
symmetry breaking and restoration. The work of [5] provides a theoretical 
framework for the study of the interplay of monopole and fermion dynamics, 
and a numerical strategy for its quantitative verification. Unfortunately, it has 
not been developed into a predictive theory, complete with equations coupling 
monopole and fermion observables and their interrelated scaling laws. Rather, 
it is just a 'physical picture' which could even be internally inconsistent. The 
speculative scenario is interesting and will be summarized here because we will 
find some support for it in our simulation data (although more conventional 
scenarios can accomodate the data as well). 

Consider massless, pointlike fermions propagating in an ensemble of U{1) 
monopoles. Suppose that the fermions experience only the U{1) gauge in- 
teraction. The fermion-monopole interaction is purely topological in origin 
because its strength is given by the product of the electric charge e and the 
magnetic charge g which are related by the Dirac quantization condition, 
eg = 2ti. This condition is renormalization group invariant, guarantees that 
the fermion-monopole interaction is always 'strong' and survives conventional 
perturbative screening, and lies at the heart of electric-magnetic duality. Sup- 
pose that the monopole network in the vacuum state of such an abelian U{1) 
model experiences a second order phase transition at a certain coupling cm 
where the monopole loops grow without bound. Let v denote the correlation 
length exponent for the transition. Suppose also that ordinary photon ex- 
change is not strong enough to cause a chiral transition near cm- Then, since 
the fermion's chirality is not a conserved property in the presence of the point- 
like monopole, the monopole transition could induce a chiral transition and 
the two point correlation function of ipip should inherit the critical exponent 
V from the monopole-loop transition. Past measurements of monopole per- 
colation have determined v with good accuracy, v ~ 0.67, which is near the 
correlation length exponent of ordinary bond percolation in four dimensions. 
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In any case, since v is distinct from the mean field value of 1/2, this mechanism 
has been cited as a potential origin of interacting fermion field theories. 

This mechanism has been explored so far in non-compact models [5-9] , where, 
however, monopole solutions cannot be found analytically, monopoles are 
(naively) expected to decouple in the continuum limit, the Dirac string costs 
action and static monopoles do not generate long-range Coulomb fields. The 
ideal testing ground for the scenario of [5] is a simple model with com- 
pact fields, a second order chiral transition and monopole percolation. When 
we begun this study the U{1) chiral transition was thought to be first or- 
derQ. These considerations lead us to consider a fermion-scalar-compact U (1) 
model, where the existence of a second order chiral transition was well es- 
tablished both analytically and numerically [14,15]. Note that this transition 
separates confinement from Higgs phases, while the QED chiral transition 
separats confining and Coulomb phases. 

From a field theory point of view, this work continues the search for an interact- 
ing field theory which is strongly coupled at short distance [16]. Considerable 
work along this line has been invested in the study of the pure scalar-gauge 
model, with compact and non-compact gauge fields. Detailed numerical stud- 
ies have supported the idea that non-compact scalar QED is trivial a la A0'^ 
by finding evidence for the tell-tale logarithms [17,18]. Analogous results were 
reported for the compact model in [19]. The simulation results support the 
longheld theoretical prejudice that there is complete charge screening and the 
renormalized gauge coupling vanishes. In addition, in non-compact scalar QED 
monopole percolation turns out to be unrelated with the Higgs transition [18]. 
These results rise the question if the chiral nature of the transition in fermionic 
models might make the continuum limit non-trivial. 

In a particle physics context, fermion-gauge-scalar models have been consid- 
ered as possible toy models for the strongly coupled standard model [20,21]. 
More recently, one such lattice model - dubbed in this context x^</'4 ^ has 
been reconsidered in [22-27] as a new mechanism for spontaneous symme- 
try breaking and dynamical mass generation. Recent results showed that the 
chiral transition line flows into the Higgs transition line, terminating with a 
tricritical point, whose critical behaviour indicates nonpertubative renormal- 
izability, consistent with a Yukawa theory. The long term aim is to find out 
whether and how the physics of the xU 04 model is different from that of the 
usual Higgs- Yukawa mechanism of fermion mass generation. For this reason, 
the triviality issue is not really central here: it might well be that the model 



^ The recent results on pure compact QED [10-12] indicate a second order transi- 
tion and suggest that compact QED could be a very interesting testing ground for 
the scenario discussed in this paper, possibly also for the full theory, if negative 7 
is chosen [13] (extended Wilson action with — 7X]p cos(20p) term). 
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provides an alternative to the usual Higgs- Yukawa mechanism, with higher 
upper bounds for its validity. 

It would be interesting if the results of this paper could be of use to further 
sharpen the understanding of the tricritical point of the xU4>4 model. However, 
this is a quenched study and its relationship with the full model can only be 
judged a posteriori: a better understanding of the xf/04 model, and/or of the 
validity of the quenched approximation, is not the primary scope of this work. 
Our aim here is to illustrate some general dynamical mechanisms possibly 
relevant to continuum physics by investigating the interplay between chiral 
phase transition and monopole percolation. We choose to do so in the simplest 
possible situation: a scalar U{1) gauge model with quenched fermions. 

This paper is organized as follows: first we review the action and phase diagram 
of the model. We continue by discussing the topological excitations of 

Abelian models, and the observables relevant for the chiral transition. Section 
4 gives an overview of the results for the topological excitations and show their 
correlations with other observables characterizing the phase diagram. These 
results confirm and extend the previous work of [28] . After this general study, 
we focus on three points along the chiral transition line: in Section 5 we present 
the results for the percolation of magnetic monopoles, in Section 6 those for 
the chiral transition. We conclude with a brief summary and discussion. 



2 Action and phase diagram 

The xU 04 model is defined by the action [22] 

S = Sy^ + Su + 
1 ^ 

^x^:^J2XxYl Vx^iUx,nXx+i^ - Ul_t,,i,Xx-^) + a-rriQY.XxXx 

^ X fj,= l X 

Su = -PY.cos{Qp) 
p 

4 

S^ = -K^Y1 {(t>lUx,n(i>x+n + h.c). 

X ^=1 

Xx are the Kogut-Susskind fermion fields, (j)^ is a scalar field with frozen length 
1 I = 1, and Ux,^ represents the compact link variable. Ux,ij,v is the plaquette 
product of the link variables Ux,^j. and r^^^; = 

The phase diagram of the quenched approximation is shown in figure 1. It has 
two phases, the Coulomb phase and the confinement- Higgs phase. The second 
is separated into two regions with different realizations of chiral symmetry. 
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Fig. 1. The phase diagram of the Higgs model (xt^04 model for amo = oo in the 
update). The chiral transition of the quenched xU4'i model is marked dotted. The 
points simulated in this work are those marked by open circles and the points N 
and E. The labels are the same as in [22,23]. 



which are seperatc phases in the chiral limit of the full model. The region 
with small k and j3 is chirally asymmetric, the one with large j3 and k is 
symmetric. These two regions are sometimes also called the confinement and 
Higgs regions, respectively, because of the limiting cases k — Q and (5 — oo. 

For /3 = the gauge and scalar fields can be integrated out exactly, leading to 
a lattice version of the NJL model [14] (point N). The point E is the critical 
end point of the Higgs phase transition line ETC, which is a tricitical point 
in the chiral limit of the full model. T is the triple point and C the phase 
transition from the confinement (at strong gauge coupling) to the Coulomb 
phase (at weak gauge coupling) in the model without the scalar field. S is the 
critical point of the spin model (/3 — oo). 

The status of the current investigation of the full model, including the depen- 
dence on the bare mass, can be found in [23-27]. Even though past lattice 

studies in the full model concentrated mostly on the tricritical point E, the 
scaling behavior along the chiral phase transition NE except E is believed 
to be NJL-like with reasonable accuracy [23]. At E the scaling behavior is 
different, e. g. the exponents seem to be nonclassical, which makes the point 
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particularly interesting. 

In this paper we study the quenched model, which corresponds to the limit 
arrio = oo in the update. The dynamics is the same as in scalar QED, or 
Abelian Higgs model, and is only considered for the measurement part. 
The qualitative features of the phase diagram are known to be the same as 
in the full model. The endpoint of the Higgs transition line E is known with 
great accuracy [19]. We expect that also in the quenched model the point 
N is described by mean field exponents, as observed in [29], probably with 
logarithmic corrections. The topological excitations of the model have been 
considered in [28]. In our study we will characterize more fully the phase 
diagram using topological excitations and we will study in quantitative detail 
the chiral transition line extending from the point N to the endpoint E, and 
its interplay with the monopole percolation. 



3 Operators 

In Abelian models with compact fields the periodic nature of the interactions 
makes it possible to rewrite the Action in terms of its topological excitations: 
magnetic monopole loops in pure compact QED, and closed vortex surfaces for 
pure scalar matter. For scalar-gauge systems there arc strings of magnetic flux, 
i. e. lattice monopoles, closed vortex surfaces and open surfaces bounded by 
monopole loops. These topological excitations, investigated in the continuum 
in [30], are amenable to numerical studies on the lattice [28,31,32]. We review 
the main definitions below. 

The study of the chiral transition will follow [33-35]. The relevant operators 
are reviewed in the last subsection. 



3. 1 Monopoles 

For the measurement of the monopole loop density we introduce the flux 
variables [31] 

Opa{.x) = e^^p + e^+p^a - Ox+a,p " O^^a £ (-47r, 47r] (1) 

where 9x,p € (— 7r,7r] is the phase angle of the link variables Ux,p and the 
physical flux 

Opa{.x) = 9pa{x) + 2Tmp„{x) e (-7r,7r] , (2) 
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where npcr{x) is an integer. If riper {x) ^ a, Dirac sheet goes through the 
plaquette. The balance of the flux entering and leaving a 3-cube defines (on 
the dual lattice) a monopole current 

2nmf,{x) = e^upaVjpa{x) . (3) 

The operator Vj, labels lattice differentiation. The vector mn{x) is the total 
flux out of the cube c at the dual site x in direction 

2nm^{x) = J2 0p= J2i^p + ^Trn^) = 2n ■ (4) 

The loops of the monopoles (on the dual lattice) are closed, because of current 
conservation 

V^m^(a;) = . (5) 

They are gauge invariant objects, whereas the Dirac sheets, apart from their 
boundaries, can be distorted. 

From this we calculate the monopole loop density, in short monopole density, 
as 




Monopole percolation is detected using an order parameter borrowed from 
standard percolation models [36]. A connected cluster of monopoles is intro- 
duced [32]: one counts the number of dual sites n which are connected with 
each other by monopole line elements, n is the size of the cluster. Note that 
this construction ignores the vector structure of the monopole currents. 

The density of the occupied bonds reads 

/y-, 1 ^max 

'Hot ^ 

where Qn is the number of clusters with size n, rimax is the size of the largest 
cluster, and ritot is the total number of connected sites. 



An order parameter for percolation is 



Its associated susceptibility is 

XM = V ({MU - (Mpe.c)') . 

We have also considered the same susceptibihty as in [32] 



Xperc 



n=4 



/rit, 



;ot 



(9) 



(10) 



3.2 Vortex sheets 



For the measurement of the vortex sheets we introduce [3,28] 
■du{x) ^(fix- Ox,u - ^x+u + aj,{x) e (-TT, tt] , 



(11) 



where ipx £ (— tt, tt] is the phase angle of 0-^ and a,j{x) is a integer. The local 
vortex sheet is now 



(12) 



Prom this we calculate the vortex density as 



(13) 



3.3 Chiral condensate 



The chiral condensate is measured in the usual way 



(14) 



with a stochastic estimator, where M is the fermion matrix. 

We also measured the logarithmic derivative of the chiral condensate [34] : 



— 



d\n{xx) 



9 In arriQ 



^ amp d{xx) 
p,^ ixx) damo 



(15) 



/3a 



This can be expressed as ratio of zero momentum meson propagators C{p — 0) 
(susceptibilities) : 



djxx) 
damn 



a(p = o), 



(3,K 



ixx) 

arriQ 



Cnip = 0) , 



(16) 



where the second is the Ward identity, which results from the chiral U{1)a 
symmetry of the staggered fermions [37]. This gives 

C^{p = 0) 



Wc note that in the quenched approximation only the connected part of Ca 
need to be considered, since fermion loops are neglected [38] . We have checked 
that our measurement of Rt^ using eq. (17) is in good agreement with the 
logarithmic derivative of the chiral condensate computed by numerical differ- 
entiation of the (xx) results. 



4 Measurement of the topological excitations - Overview 



Fig. 2 shows a overview of the results on the 6^ lattice. For /? = 0.00 (a) 
(NJL hue) all observables are very smooth, only the ratio Mperc and the chiral 
condensate show a somewhat steeper descent. Monopole and vortex density 
decrease slowly for increasing f3. More to the right on the NE line for f3 = 
0.60 (b) a similar behavior is observed but the transition region is somewhat 
smaller. 

For /? = 0.90 (c) (crosswise to the ET line first order) a steep descent in all 
observables at /t = 0.24(1) can be observed, as this is typical for first order 
phase transitions. Obviously the monopoles and the vortices are sensitive to 
this phase transition. 

At (3 = 1.50 (d) (crosswise to the TS line first order) a phase transition 
from Coulomb to Higgs phase can be observed. At this phase transition only 
the vortex density and the link energy are sensitive, which vary rapidly at 
K, — 0.17. The other observables stay nearly constant. 

The phase transition CT from confinement to Coulomb phase is shown for 
K = 0.05 (e). At the phase transition at (3 = 1.01(1) the monopole density, 
Pm and the chiral condensate steeply decrease. Link energy and vortex density 
are insensitive to this phase transition. 
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Fig. 2. Overview of the topological observables pm, Mpcrc and pv together with El 
and (xx) on the 6^ lattice at the lines NE (a+b), ET (c), TS (d) and CT (e). The 
errors (not shown) are nearly always smaller than the symbols. 
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Fig. 3. Monopole density pm and vortex density pv for (a) /3 = 0.848 ~ /3e and (b) 
P = 0.30. 

So, we have the following picture. Consider the three regions of the phase 
diagram: Coulomb phase, Higgs region and confinement region. In the con- 
finement region the monopole and vortex density are large and in the Higgs 
region both are small. In the Coulomb phase the vortex density is large and 
the monopole density small. At the first order phase transtions the observables 
which are large in one phase and small in the other show a steep decrease, 
whereas the others are insensitive. 

The order parameter of percolation Mperc follows (essentially) the monopole 
density but with more rapid changes. This could be understood in analogy to 
bond percolation. If the monopole density is high, percolating clusters are very 
likely, whereas for small monopole densities they are very unlikely. The first 
percolating clusters occur at some intermediate coupling. For large /3, along 
the line of first order phase transitions, this seems to coincide with the phase 
transition of the monopole density. Such interplay of monopole condensation 
and percolation has been investigated before in pure compact QED in [39] . 

At the fine NE, which is second order in the full theory, the situation is less 
obvious. 

At the endpoint E of the Higgs phase transition, (3 — 0.848, and pv 
show a steep descent, which gets increasingly steeper for increasing lattice size 
(Fig. 3a). The derivative of these observables show an increasing maximum. 
We have not done a scaling investigation but the data suggest a divergence. 

For P — 0.30 (Fig. 3b), still along the line NE, both densities show no max- 
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K 



Fig. 4. Percolation susceptibilities Xpcrc (solid) and xm/10 (open symbols) for 
(3 = 0.30. The curves are interpolations made with the use of reweighting tech- 
niques. 

imum in the derivative. Moreover, no volume dependence was observed. Ob- 
viously they do not get critical in the region where the chiral condensate has 
its largest variation. 

For P — 0.00 the picture is very similar to that at /3 = 0.30. 

The three /3 values 0.00, 0.30, and 0.848 have been selected for the systematic 
investigation which we describe in the rest of this paper. 



5 Magnetic Monopoles in the vicinity of the NE line 

In this Section we measure the critical line for the percolation of the magnetic 
monopole clusters, and its critical exponents. An accurate measurement of the 
critical line is obtained thanks to the combined use of reweighting techniques, 
and the two definitions of monopole percolation Xperc and xm- 

The quahty of our measurements is demonstrated in Fig. 4 where we contrast 
Xperc and Xm at /3 = 0.30 on lattices of different size. Both susceptibilities 
show a clear maximum for all /3's in the vicinity of the line NE. The peaks are 
very close, and the volume corrections have different signs. 

We have determined the precise position of the maximum of Xperc and xm 
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Table 1 

Position and height of the maximum of the percolation susceptibility for different 
(3 and lattice sizes along the line NE. The values are given for both definitions xm 
and Xperc- For the infinite volume extrapolation of Kmax = 1-3469, k, = 0.8875, 
and K = 0.2714, respectively) also the density of the occupied bonds p„i and the 
monopole density pm is printed. For comparison, our best estimate for the chiral 
transition in the infinite volume is given in the second column (see section 6). 





chiral, c 


lattice 




Xm{L) 




Xperc (i) 


Pm 


Pm 


0.00 


1.23(3) 


6^ 
84 
12^ 
16* 

20-* 


1.307(4) 

1.3179(6) 

1.3303(5) 

1.3358(6) 

1.3395(3) 


156.(3) 
342.(2) 
941.(9) 
1865.(13) 
3054.(31) 


1.363(5) 

1.3525(6) 

1.3474(4) 

1.347(1) 

1.3467(4) 


16.11(5) 
32.00(4) 
80.9(4) 
154.(7) 
253.(3) 


0.1132(1) 

0.11323(2) 

0.11322(1) 

0.11322(1) 

0.11321(1) 


0.5283(6) 

0.52852(9) 

0.52848(5) 

0.52847(7) 

0.52841(4) 


0.30 


0.87(2) 


64 
8* 
12* 
16-* 


0.861(4) 
0.8671(5) 
0.8756(4) 
0.8812(4) 


165.(3) 
360.(2) 
989.(5) 
1937.(13) 


0.903(3) 
0.8926(8) 
0.8886(8) 
0.8882(5) 


15.6(2) 
31.5(1) 
80.3(4) 
153.7(7) 


0.1074(3) 
0.10746(3) 
0.10749(2) 
0.10750(1) 


0.494(1) 
0.4946(2) 
0.49474(9) 
0.49477(7) 


0.848 


0.26:0.29 


84 
12^ 
16* 
20* 


0.2678(1) 
0.2692(1) 
0.2701(1) 
0.27065(5) 


574.(4) 
1437.(15) 
2758.(34) 
4501.(53) 


0.2732(1) 
0.27204(7) 
0.27191(8) 
0.27168(4) 


22.35(6) 
59.5(3) 
120.(1) 

206.1(8) 


0.0727(2) 
0.07257(5) 
0.07259(7) 
0.7257(2) 


0.3180(7) 
0.3174(2) 
0.3175(3) 
0.3174(1) 



with the multi-histogram method. The numerical results for Xperc and xm, 
together with other magnetic observables, are summarized in table 1. In Fig. 5 
we show the position of the peak of the two susceptibilities as a function of 
1/L for /5 = 0.00. We superimpose a power law fit with different amplitudes 
and exponents and a common critical point. Again we see that, by increasing 
volume, Xperc approaches the critical point from above, xm from below. They 
provide upper and lower bounds, hence an accurate determination, of the 
critical coupling for the percolation of magnetic monopole clusters. 

A determination of the critical exponent of comparable accuracy requires a 
careful consideration of scaling violating effects, and will be presented else- 
where. For our present purposes a precision comparable to that we achieved 
at the chiral transition would suffice (5-10 %). 

Fig. 6 shows the scaling behavior of the maximum of Xperc with the lattice 
size. Note the linear rise with size in this double logarithmic plot, supporting 
power law scaling. Only the data on the 6^ lattice displays small deviations. 
In a fit with the finite size scaling equation 

Xperc,max(i^) OC L^p-c/^p- (ig) 

we measure for (5 — 0.00 and (5 — 0.30 the nearly identical value of 7perc/ ^'perc — 
2.28(2). The subscript 'perc' is added to indicate the definition of the ex- 
ponents at the percolation transition. For comparison, the measurements in 
quenched non-compact [6] and in full non-compact QED [5] yields quite sim- 
ilar values of 7perc/i^perc = 2.24(2) and 7perc/i^perc = 2.25(3), respectively. For 
/3 = 0.848 one obtains a somewhat larger value of 7perc/t'perc = 2.42(1). 
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Fig. 5. Position of the peak of the percolation susceptibihty Xperc (circles) and xm 
(triangles) for P = 0.00. 




fit: A^L^/" P"™: 
O:/? = 0.00 

y/u^,,,= 2.269(6) 
A= 0.286(4) 
X^do,= l. 985187 

^:(^ = 0.30 

y/v^,,= 2.288(8) 
A= 0.271(5) 
X^do,= 1.947799 

□ :/? = 0.848 
r/i/pe.c= 2.425(5) 
A= 0.1443(18) 
X%or=0. 369933 



Fig. 6. Scaling behavior of the maximum of the percolation susceptibility Xperc for 
(3 = 0.00, (3 = 0.30 and (5 = 0.848. In the fit only data with L > 8 have been 
considered. 
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The maximum of xm shows larger deviations from the leading finite size scaling 
behaviour for small lattices. If we use only lattices with L > 12, we get for 
/3 = 0.00 and (3 = 0.30 an exponent 'Jm/i'm = 2.32(3) and for (3 = 0.848 
Im/i^m = 2.23(4). The subscript M indicates the exponents which have been 
determined using xm- 

The 7/z/'s from the two definitions are separated by less than two standard 
deviations at /5 = 0.00 and /3 = 0.30. At (3 = .848 the disagreement increases. 
Larger lattices, and/or a careful consideration of correction to scaling terms 
should bring the results from the two definitions closer. Within this uncertainty 
all numbers of of the percolation transition are compatible with 2.3(1). 

We can contrast these results with those of pure random bond percolation 
[40]. For pure random bond percolation the critical density is Pc — 0.161. In 
our model the density Pm,c at the percolation phase transition depends on /3, 
and it is generally lower. Our result is close to that of site bond percolation 
7/1/ = 2.094(2) [41,40]. Considering the vector character of the monopole 
current, we should not expect that the two models are in the same universality 
class, as pointed out in [6]). 

To determine /Spcrc/ ^'pcrc we investigated the scaling with volume of M^crc at the 
(infinite volume) critical coupling. The results are very sensitive to the precise 
value of Kc- We estimated anyway /Sperc/^'perc = 0.85(3) in good agreement 
with the scaling relation 



The result for quenched non-compact QED is /3perc/i^perc = 0.88(2). 

The exponent /3perc was determined from Mperc in the phase with percolation, 
using the results not distorted by finite size effects. The result for f3 = 0.00 is 
shown in figure 7. Here again we fixed the critical coupling. The fits give /9pcrc = 
0.52(2), where the error was determined by varying the critical coupling in the 
interval k,c — 1.346 . . . 1.348. We repeated this for /3 = 0.30 and obtained a 
Pperc — 0.50(4), smaller but compatible within errors. (/3 = 0.58(2) in non- 
compact quenched QED [6].) 

Together with our determination of /Jperc/'^perc = 0.85(3) we get the expo- 
nent for the correlation length Uperc — 0.61(4) to be contrasted with that 
of non-compact QED i/perc = 0.66(3), and of pure random site percolation 
^'pcrc = 0.69(1). A more careful analysis, considering corrections to scaling, 
would be necessary to assess the significance of these differences and to decide 
if percolation is in the same universality class for all /3 < 0.848. Differences, if 
any, are anyway small. 




(19) 
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1.25 1.30 1.35 



Fig. 7. Scaling behavior of Mperc for (3 = 0.00. In the fit only data with L > 16 in 
the interval k = 1.30. . . 1.33 have been considered. The critical coupling was fixed 
to Kc = 1.3469. 

Note that for [3-^ = 0.848 the percolation transition at k = 0.2714(3) is distinct 
from the Higgs phase transitions which occurs at Kc = 0.26333(1) [19,25,26]. 
In Fig. 8, we see that all the observables are very sensitive to the Higgs phase 
transition (note the sharp peak of the susceptibility of the link energy), with 
the exception of the susceptibility Xpcrc of the percolation. This shows that at 
the endpoint, the Higgs phase transition and the percolation of the monopoles 
is uncorrelated. So the percolation of the monopoles probably can not help us 
understand the tricritical point in the dynamical model. 



6 The chiral transition (NE) line 

We consider first the behavior of the susceptibihty ratio. The discussion here 
aims at illustrating the general trend along the critical line. In the second 
subsection we cross check the results by fitting the data to the equations of 
state. That includes a more detailed discussions on the selection of the scaling 
window, and a careful error analysis. 
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Fig. 8. Susceptibility of the monopole percolation Xperc and the link energy xe^, 
as well as the order parameter of the percolation Mperc, the monopole density pm 
and the chiral condensate (xx)('^"^o = 0.04) as function of «; on a 12^ lattice at 
(3 = 0.848 at the endpoint of the Higgs phase transitions. 



We measure the ratio of R^^ as a function of arriQ: at the chiral transition i?^ 
should be independent on aruQ if the transition is described by a simple scaling 
law. Logarithmic corrections modify this behaviour. 

For f3 = 0.00 (Fig. 9a) our data indicate a chiral transition at k^. = 1.23(3) 
with = 0.35(5). This value compares well with the expected mean field 
result Rt^^Hc) = 1/5 = 1/3. Logarithmic corrections, if any, are small, as we 
will see later on. 

For 13 = 0.30 (Fig. 9b) the results suggest = 0.87(2) with i?^(/€c) = 0.45(10). 
This corresponds to 6 = 2.2(6). Logarithmic corrections a la Nambu-Jona 
Lasinio would predict 

i?,(«e) = l/(3 + l/log((xx))) (20) 



resulting in a slight upward trend of the critical ratio, still compatible with the 
data. We comment more on this point in the next subsection. The curves for 
aruQ > .02 — .04 starts bending downwards, after the initial rise: masses larger 
than .04 cannot be used safely to assess the critical behaviour. In the next 
subsection we discuss in detail how to select the appropriate scaling window. 
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Fig. 9. as function of amo at (a) (3 = 0.00, (b) /? = 0.30 and (c) /? = 0.848 ~ 
(Solid, dashed, dotted lines) are for (16^, 12^, 8*^) lattices. Approximate values for 
the critical point are indicated with horizontal lines. 
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/? = 0.848 

12": 

□ :/c=0.28 
*:/f = 0.27 
n:/c=0.265 
A:/c=0.26 
o:/c=0.255 

16": 

■:/f = 0.28 
+ :/c = 0.27 
*:/f=0,26 





.01 .02 .03 .04 .05 .06 

Fig. 10. Behavior of the quadratic tt meson mass as function of amo. In the broken 

phase a straight hnc throTigh the origin is expected. The hncs indicate a fit to 
the 16^ data. For k = 0.27 we observed large finite size effects, which make the 
extrapolation unreliable. 

For (3 = 0.848 (point E) the Higgs phase transition is at Kc — 0.263. We have 
shown above that the Higgs transition is separate from monopole percolation, 
and here (Fig. 9c) we see that at this k a conventional second order chiral 
transition is incompatible with the data. The data in Fig. 9c would rather 
suggest a chiral transition in the interval [.28 - .30], even on the right of the 
percolation transition. But this behaviour is far from clear. 

As an alternative strategy to locate the chiral transition a,t P — 0.848 we 
investigated the scaling behavior of (am^)^ as a function of amo (Fig- 10). 
For K = 0.255 and k, = 0.26 the tt meson mass scales corresponding to the 
PCAC relation {aniT^Y oc airiQ. For larger k, especally for those k between the 
transitions, deviations from the PCAC relation can be observed. The data for 
K = 0.27 on the 12^ lattice might suggest a straight hne joining the origin and 
the two smaller points. But the data for amo = 0.01 on the 16^ lattice indicate 
that we have to be aware of finite size effects and the data on the larger lattice 
favor a violated PCAC relation. These ambiguities in the extrapolations are 
not uncommon in the numerical studies of phase transition. We conclude that 
also the PCAC analysis does not reliably locate Kc of the chiral transition at 
(5 = 0.848. 

Because of the unusual behavior of the chiral condenstate the equation of state 
also fails. So we cannot decide from this data if the Higgs and chiral transition 
meet at the endpoint, as observed in the dynamical theory [23,26]. 
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6.2 Equation of State 



The results of the susceptibihty ratio pose at least two questions: 1) is the 
critical behaviour ai (3 — 0.30 really different from the one at /3 = 0.00? 2) if 

yes, and the data at f] = 0.30 are described by a power law scaling, what are 
the other critical exponents? - in particular, we are interested in the relative 
position of percolation and chiral transition, and in their correlation length 
exponents u. 

The standard way to address these questions is to exploit the chiral equations 
of state (EOS). First we briefly discuss the results at /3 = 0.00, which we 
shall use as a term of comparison, and then we will devote more time to the 
discussion at P — 0.30. 

The EOS framework is borrowed from ancient studies of ferromagnetic tran- 
sitions, and, as such, uses data in the thermodynamic limit. Happily, this 
approach requires that we simulate the model within its scaling window but 
not necessarily directly at the critical point. Consequently it is possible to 
work with data free from finite volume effects (which won't be possible at the 
critical point where a finite size scaling analysis would be mandatory). Both 
the selection of the scaling window, and the control over finite size effects are 
discussed below. 

Recall the possible critical behaviour, and the relative equations of state. If 
scaling holds, the data are described by a 'standard' equation of state 

amo^{xx)'fit/{xxy^^-), (21) 

where t is the reduced coupling (in our case t — k — Kc) and whose first order 
approximation - used in actual fits - reads 

amo = a^(«; - «;e)(xx)'-'/''^ + b"" {xx)' ■ (22) 

If the theory is trivial, the most natural candidate to describe the data is an 
equation of state a la Nambu-Jona-Lasinio. We will use the following form, 
which is motivated by the leading term of the 1/N expansion [42]: 

amo = a^^^(« - k,){xx) + ^^^'^(xx)'(log(xx)/«) (23) 

and we will just comment on its possible generalizations [46,47] (in practical 
analysis log(xx)/s will be replaced by log(xx) + c). 

In many cases we used two different fitting procedures. 
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Fig. 11. Results of the power law fits described in the text at /3 = 0.00. The arrows 
indicate the fitted interval. 

For the first one we numerically invert eq. (22) for (xx) then make a 
minimum fit for the measured chiral condensate, taking into account the 
errors. We label the results obtained with this procedure [Fl]. 

The second one [F2] is a least squares procedure which minimizes (amo — 
am((xx)))^, am{{xx))T being the r.h.s. of either eq. (22) or eq. (23). We 
then compTitc the errors by jack-knifing the results obtained by discarding 
one point at a time. The quality of the fit is estimated hj Q = (mo — 
m((xx)))^/(A(m((xx))^)A^points)- For a good fit we expect Q close to one. 



6.2.1 (3 = 0.00 

As mentioned above, we expect that at (3 = the model reduces to a lattice 
NJL model. Recent numerical studies of four dimensional NJL models include 
[42-44]. 

At (3 = 0.00 we restricted ourselves to k, = (1.15, 1.20, 1.25). A power law fit 
with 5 free parameters which uses the data at bare masses .02 and .01 from 
the 16^ lattice, the data at .015 from the 12^ lattice is shown in Figure 11. 
The quality of the fit [Fl] is nice, the exponents S = 3.06(19), = 0.56(4) 
are all consistent with mean field theory. The critical couphng Kc is 1.228(7). 

We have then included logarithmic corrections, initially without scale (the 
parameter s = 1 in the EOS). We used data from the 12^ lattice, and the 
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procedure [F2]. In the same interval we obtain Kc = 1.259(4). Adding one 
further point at k = 1.30. arriQ = .02 moves the central values to Kc = 1.258. 
The inclusion of a free scale parameter in the fits moves the critical point to 
Kc — 1.248, closer to the power law results. We have checked that we can 
discard the points with bare mass .02 without altering the results. We can 
also constrain the exponents to their mean field values without appreciably 
degrading the quality of the fits - that corresponds to 'cancelling' the logs in 
the NJL equation of state. Kc is then 1.224(8). 

In conclusion, the data at P — 0.00 are well described, as expected, by a mean 
field critical behaviour, with Kc in the range 1.22 — 1.26, in agreement with 
the ratio analysis. This range of critical couplings is well separated from that 
of monopole percolation as can be seen from Table 1 and Fig. 5. Logarithmic 
corrections, if any, are small, their inclusion improves the quality of the fits, 
but does not alter Kc nor do they allow fits on a wider interval. These critical 
couplings are significantly different from those of the percolation transition. 



6.2.2 (5 = 0.30 

Our first task is to select a reliable data sample for our fits. Prom the results 
for we know that masses larger than .02 are too heavy. On the other 
hand, small masses, which require large lattices, are only available for three 
k's. All in all, our candidate window for fits narrows down to the nine points 
at K — .85, .87, .90 - we occasionally checked the stabihty of our results by 
including the two arriQ — .02 points at k = .80 and .95 but we did not venture 
any further. By contrasting data in this range obtained on a S*^, 12"^, and a 
16"^ lattice, wc checked that finite volume effects are already small on a 8^ 
lattice, and that the difference between results on a 12^ and on a IQ* lattice is 
statistically not significant (the worst case is amo = 0.005 were the data are 
separated by two standard deviations). We used the data on the 16^ lattice 
for our fits. 

The results from [Fl] applied to the power law form are shown in Figure 12. 
The stability of these results was checked by performing the [F2] fits on the 
same range of parameters. The results arc shown in Table 2, first line. We get 
Kc = .874(6), S = 2.12(15), and = 0.88(5). (These exponents are very close 
to those of quenched, non-compact QED: S = 2.12(1), f3^ = 0.86(3), 7 = 1.0 
[45,35].) Assuming scaling, these critical exponents give u = .695. We have 
verified that the fit is stable against the inclusion of the two points at arrio — 
.02, K — .80 and .95. Finally, since we want to contrast the position of the 
chiral transition with the position of the peaks of the monopole susceptibility, 
we have also tried fits constraining the critical k = .882, at the the lower 
bound of Kc percolation. The results for the critical exponents (second line of 
Table 2) are compatible with the previous ones and give u — .73. The quahty 



22 



A 
X 
X 
V 




/? = 0.30 

O:amo=.02 (16") 
□ :amo=.01 (16*) 
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X%of=2.018288 



Fig. 12. As Figure 11, but /3 = 0.30. This fit correspond to the first hue in table 2. 



Table 2 

Results of the power law EOS (eq. (22)) fits at /3 = 0.30. 









6 




Q 


2.54(4) 


1.58(36) 


.874(6) 


2.12(15) 


.999(68) 


1.068 


2.19(28) 


1.22(7) 


.882 [Fixed] 


1.94(2) 


.932(47) 


1.79 


2.00(40) 


1.08(8) 


.8875 [Fixed] 


1.83(3) 


.893(82) 


2.96 



parameter Q (1.79) is still acceptable. Constraining Kc deeper into the critical 
interval spanned by the critical k for percolation further degrades the quality 
of the fit, which remains neverthless tolerable. We have also searched for the 
best exponents assessing by eye the quality of the scaling plots built following 
eq. (21). The best plot picks Kc = -868, 6 = 2.26, = .8, giving u = .65. 

We have also tried to fit the data with a logarithmic improved mean field 
equation. The results of the fits are displayed in Table 3. (Note, that this fit 
has less free parameters.) We see that the only fit of quality comparable to 
the power law fits requires a variable scale. In this case, the critical coupling 
is smaller than that from the power law fits: k,c — .866(2) to be contrasted 
with Kc — .874 (6). The shift in the critical coupling with respect to the 
power law results can also be appreciated by performing a constrained fit 
with Kc = .882, which fails completely (third line of Table 3, Q = 83.26): this 
definitively shows that a NJL critical behaviour cannot be associated with the 
percolation of magnetic monopoles. 
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Table 3 

Results of the logarithmic EOS (eq. (23)) fits at /? = 0.30 



^NJL 


^NJL 




c 


Q 


2.46(6) 


-3.68(16) 


.861(1) 


[Fixed] 


6.94 


2.54(5) 


-6.12(84) 


.866(2) 


.69(14) 


1.18 


2.81(19) 


-13.88(2.01) 


.882 [Fixed] 


1.23(13) 


83.26 



.70 
.65 
.60 
.55 
.50 
.45 
.40 
.35 
.30 
.25 
.20 



: /? = 0.30 

-- a:K = 0.9 
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Fig. 13. in the critical region. Dotted is for the 16^ lattice, dot-dash for the 12^. 
The dashed line is i?jr(/^c) = 1/S from the power law fit Table 2, the solid line is 
Rk{kc) from the logarithmic fit, first line Table 3. 

We can also study the sensitivity of the position of the critical point to the 
equation of state by inspecting Rt^. This gives a clearer indication of the role 
played by the precise location of the critical point in determining the critical 
scaling. In Fig. 13 we show the results for in the critical region constrasted 
with 1/S (sohd Hne) from our best power law fit, and with Rt:{k.c) from the 
logarithmic fit without scale. We see that the i^^'s corresponding to the two 
hypotheses (logarithms and power law) fall in the k, interval predicted by the 
fits - between .85 and .87 for log fits — .86), and between .87 and .90 for 
power law — 



Summarizing: If the correct results (i.e. those describing the physics of 
the chiral limit on infinite volumes) are given by the logarithmic fits, second 
line Table 3, the percolation transition and the chiral transition are distinct. 
In this case we have not learned anything about a possible role of magnetic 
monopoles in the existence of strongly coupled theories. Simply, the chiral 
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transition woTild be trivial and magnetic monopoles would have nothing to do 
with it, as it is for j3 = 0.00 and probably also for /3 — 0.848. 

If the correct results are given by power law fits, Table 2, the chiral transition 
coincides within errors with the percolation transition. The critical exponent 
V for the chiral transition is in the range .65 - .73, to be contrasted with 
i^perc = -61(4), as we determined it in section 5. Taking into account the 
possibility of further corrections, the exponents v for the two transitions are 
compatible and ~ 0.65. 

Note again these results differ from the ones in the full model where the crit- 
ical behavior along the line NE seems the same as the one of the trivial NJL 
point. It is not in the scope of the present paper to investigate the relation- 
ship between the quenched and full model, but this is certainly an interesting 
subject for further investigations. 

Finally, we might wonder about the sensitivity of the quality of the logarithmic 
fits and/or the value of the critical coupling to further corrections. Typical 
corrections include powers of the logarithms. Although such corrections are 
not there, up to 0(1/A^^), they cannot be excluded in general. Our model, 
however, could be in the Yukawa universality class [46,47]. In this case, one 
can conceive of sizeable variation, even in sign, of the exponents and the 
simplest parametrization of such crossover behavior would suggest that the 
powers in the logs should be kept free. It cannot be excluded that more general 
logarithmic equations of state fit the data well, and certainly this would make 
the study of the interrelation between the chiral transition and monopole 
percolation more subtle. 



7 Summary 

We have investigated the topological excitations in the phase diagram of scalar 
QED. We confirmed with great precision the picture suggested by [28] and see 
clear signals for first order phase transitions in the topological observables. 

We have investigated in detail chiral transition and monopole percolation for 
three points: 

• 13 — 0.00: The chiral transition is described by mean field exponents, pos- 
sibly with small (logarthmic) corrections. The percolation transition is well 

separated from the chiral transition. 

• (3 = 0.30: Two possibilities are consistent with the data: 1*^* a situation 
analogous to that at /5 = 0.00, but with large logarithmic corrections; 2°*^ 
percolation and chiral transition coincident, with the same critical exponent 
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V. 

• j3 = 0.848: The only clear statement here is that the Higgs transition and 
monopole percolation are distinct. We have not been able to measure with 
confidence the chiral transition, which has proven hard to study with con- 
ventional numerical methods. 

Since the scalar field increases the order in the gauge fields and favours chiral 
symmetry, separate chiral and monopole percolation transitions are not unex- 
pected: it might well be that monopoles percolate in the symmetric phase, but 
their tendency to break chiral symmetry is overcome by the ordering effects of 
the scalar fields. This seems to be the case ai (3 = 0.00, and it is also possible 
at /9 = 0.30 (if the first possibility is realized). In this case chiral transition is 
most probably mean field like along the whole line NE (without the point E), 
as it seems to be in the full model. Moreover, the chiral transition is separated 
from the percolation transition: when this is the case monopoles should be 
irrelevant to the dynamics of the chiral transition in the large volume, con- 
tinuum limit. Except for the point E (we reiterate that we don't have any 
conclusion on the relative position of monopoles and chiral transition at the 
point E) this goes hand in hand with the transition being trivial. Still, the be- 
haviour at /3 = 0.30 remains different from that at /3 = 0.00, and and it would 
be interesting to find a coherent scenario accomodating these observations. 

The scenario of [5] requires instead a chiral transition coincident with the per- 
colation transition, and sharing the same correlation length exponent. When 
this occurs, the dynamics of the chiral transition should inherit the character- 
istics of that of the magnetic monopoles, which survive the non-trivial con- 
tinuum limit. This could be the case at /3 = 0.30. The scenario of [5] receives 
then some support from our investigation. 

Important questions concern the relevance of these results for the full model, 
and, in general, the sensitivity to the number of flavours. If the scenario of [5] 
is realized, the situation is different from the one observed in the full model: If 
this is the case, fermion screening plays an important role in this model as it 
presumably does in the gauged Nambu Jona Lasinio model: there the ladder 
(quenched) approximation predicts a non-mean field scenario [48,49], while 
the full model is presumably trivial [50]. These observations do not detract 
from the illustrative value of our results at /5 = 0.30, but, of course, prevent 
us from extending any of our conclusions to the full model at this stage of 
our investigation. If the full and quenched models really have different critical 
behavior along the fine NE, then we are challenged to understand something 
unexpected and special about fermionic screening in this model. In the present 
quenched model, the dynamical scalar field in the configurations completely 
screens the electric charge, producing a vanishing renormalized gauge coupling 
in this sector of the model. If the fermions were also treated dynamically, they 
would add nothing qualitatively new in terms of charge screening. Therefore, 
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the fermion-monopole interaction would have to have some essential ingredient 
which renders the monopoles ineffective in driving chiral symmetry breaking 
in the full model. Simulations of non-compact QED with the number of flavors 
Nf varying from 2 to 32, do not show such an effect — there chiral symmetry 
breaking and monopole percolation were coincident, within ample errorbars, 
for all Nf [5]. Better simulations of the full and quenched models with smaller 
bare fermion masses and larger lattices are needed now. 

We should define better numerical strategies so to show more clearly the corre- 
lation, or lack thereof, between chiral and monopole observables. It is however 
evident from our data, that for (3 = 0.00, when the chiral and monopole 
transitions are clearly distinct, the chiral transition is clearly trivial and is 
not correlated with the percolation of monopoles. When the transitions grow 
closer for f3 = 0.30, so do the correlation length exponents u (or the effective 
critical exponent u in the case of logarithmic triviality) . The transitions may 
be coincident, with the same critical exponent i'. 
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